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respectively,  indicates  that  almost  uniformly  k  *  k  ,  the  latter 
values  being  typically  ca  10-lCr  fold  smaller.  Tnese  rate  discrepancies 
are  chiefly  reflected  in  the  activation  entropies,  although  partly 
compensated  by  differences  between  the  experimental  and  calcualted 
activaton  enthalpies.  Although  the  values  of  k  /k  orr  depend  somewhat 
upon  the  reaction  environment  as  determined  by  Efienafure  of  the 
coordinated  ligands  and  the  metal  surface,  they  are  approximately 
independent  of  the  magnitude  of  k  or  the  dirivng  force.  Taken 
together,  these  findings  support  £Be  notion  that  the  major  origin  of  the 
discrepancies  betweek  k  and  ^corr  is  associated  with  changes  in  local 
solvent  structure  when  forming  theprecursor  state.  Nonadiabaticity  may 
contribute  importantly  by  necessitating  that  the  reacting  centers  be  in 
very  close  proximity.  Reactions  at  "hydrophilic"  metal  surfaces,  such 
as  lead  and  gallium,  that  are  known  to  strongly  adsorb  water  molecules 
yield  remarkably  similar  theory-experiment  disparities  to  those  seen 
with  cationic  coreactants  in  homogeneous  solution.  The  markedly  closer 
agreement  obtained  for  some  reactions  at  mercury  is  attributed  to  the 
relatively  mild  perturbation  exerted  by  this  metal  surface  upon  the 
local  solvent  structure. 
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THE  DRIVING-FORCE  DEPENDENCE  OF  ELECTROCHEMICAL 
RATE  PARAMETERS:  ORIGINS  OF  ANODIC-CATHODIC 
ASYMMETRIES  FOR  METAL  AQUO  REDOX  COUPLES 

Joseph  T.  Hupp  and  Michael  J.  Weaver * 

Department  of  Chemistry,  Purdue  University 
West  Lafayette,  Indiana  47907 

Abstract 

^  The  consequences  of  differences  in  the  intramolecular  force  constants 

and  the  ionic  entropies  between  the  oxidized  and  reduced  states  of  aquo 

redox  couples  upon  their  electrochemical  kinetics  are  examined  as  a  function 

of  the  driving  force.  A  generalized  harmonic  oscillator  model  is  utilized 

that  involves  estimating  the  activation  barrier  from  the  individual  force 

constants  in  the  oxidized  and  reduced  states  rather  than  employing  average 

("reduced")  values.  Noticeable  asymmetry  in  the  anodic  and  cathodic  Tafel 

plots  is  predicted  for  redox  couples  having  large  (ca.  twofold)  differences 

in  force  constants,  the  plots  being  markedly  more  curved  at  anodic  over- 

potentials >'^he  calculated  plots  are  in  reasonable  agreement  with  experimental 

Tafel  plots  determined  previously  for  Cr(0H_)  3+//2+  and  Eu(0H_)  3+/,2+  at 

<£t>  z  n 

the  mercury-aqueous  interface.  Comparisons  are  also  presented  between 

3+/2+ 

electrochemical  "ideal"  activation  parameters  determined  for  Cr(0Ho), 

as  a  function  of  anodic  and  cathodic  overpotentials  and  the  corresponding 

parameters  calculated  from  structural  and  thermodynamic  data.  Substantially 

2+ 

smaller  activation  enthalpies  are  observed  for  Cr(0Ho),  oxidation  than 

3+  2  “ 

for  CrCO^)^  reduction;  both  are  close  to  the  theoretical  predictions. 

The  anodic  activation  enthalpies  approach  zero  at  moderate  overpotentials 

even  though  the  activation  entropies  remain  large  and  negative.  Parallels 

are  also  drawn  with  corresponding  results  and  data  interpretations  of  the 

driving-force  dependencies  for  related  homogeneous  reactions. 


Introduction 


We  have  recently  reported  and  discussed  an  anomalous  driving-force 

dependence  of  rate  constants  for  the  outer-sphere  electrochemical  reduction 

and  oxidation  involving  metal  (III)/(I1)  aquo  couples  of  chromium,  vanadium 

1  2 

and  europium  at  the  mercury-aqueous  interface.  ’  The  metal (III)  reductions 

exhibit  linear  variations  of  the  logarithm  of  the  work-corrected  rate 

constant,  ^corr>  with  electrode  potential  (i.e.  Tafel  plots)  over  the  entire 

accessible  range  of  overpotentials  (ca.  300  to  600  mV,  depending  on  the 
1  3 

reaction).  ’  In  contrast,  the  work-corrected  Tafel  plots  for  the  corres¬ 
ponding  metal (II)  oxidations  in  hexafluorophosphate  or  perchlorate  electrolytes 

are  sharply  curved  at  moderate  overpotentials,  vielding  variations  of  loc  k 

6  corr 

with  electrode  potential  at  higher  oxidation  driving  forces  that  are  unexpectedly 

small1  in  comparison  with  the  predictions  of  the  usual  harmonic  oscillator 
4 

treatment.  Similarly  weak  dependencies  of  the  rate  constants  upon  the 

thermodynamic  driving  force  have  been  noted  for  a  number  of  homogeneous 

2  5  d 

bimolecular  reactions  involving  oxidation  of  M(II)  aqua  complexes. 

The  observed  marked  asymmetry  of  the  anodic  and  cathodic  Tafel  plots 
contrasts  the  essentially  symmetrical  behavior  predicted  by  the  conventional 
theoretical  relationships.1  We  have  speculated  that  the  major  source  of 
this  asymmetry  is  associated  with  major  differences  in  the  short-range 
solvent  polarization  between  the  oxidized  and  reduced  aquo  species  ,  possibly 
associated  with  a  rearrangement  step  either  before  or  after  electron  transfer.1 
However,  the  recent  acquisition  of  reliable  Raman  spectral  measurements  for 

9 

both  M(III)  and  M(II)  aquo  complexes  along  with  the  emergence  of  new 

metal-ligand  bond  length  data  1^  has  prompted  us  to  reexamine  this  problem. 

Specifically,  these  data  indicate  that  the  force  constants  for  the  totally 

symmetric  metal-aquo  stretching  mode  are  expected  to  differ  substantially 

(ca.  twofold)  between  the  oxidized  and  reduced  forms  for  couples  such  as 

o 

Cr  which  feature  large  changes  ('0.15  to  0.2  A)  in  the  metal-ligand 

bend  lengths,  _a .  Such  force  constant  differences  will  yield  noticeably 


different  free-energy  parabolas  for  the  oxidized  and  reduced  forms  for  these 
reactions  since  a  large  portion  of  the  free-energy  barrier  is  associated 
with  inner-shell  (i.e.  metal-ligand)  reorganization.  Such  asymmetry  in  the 
free-energy  curves  for  individual  redox  couples  will  often  cancel  for 
homogeneous  reactions,  especially  those  involving  a  structurally  similar 
pair  of  redox  couples.  However,  it  can  in  principle  yield  significantly 
asymmetric  activation  energy-driving  force  [and  hence  log  ^corr“  potential] 

4 

plots  for  electrochemical  reactions  at  anodic  and  cathodic  overpotentials. 
The  extent  of  such  asymmetry  has  been  shown  to  be  only  minor  at  small 
overpotentials.^  We  have  stated  previously  that  this  factor  is  much  too 
small  to  account  for  the  observed  Tafel  plot  asymmetry.^-  However,  more 
detailed  consideration  of  the  inner-shell  reaction  energetics 
on  the  basis  of  the  recent  structural  data  leads  us  to  revise  this 
conclusion. 

The  purpose  of  this  communication  is  to  examine  in  quantitative 

terms  the  predicted  influence  of  differing  metal-ligand  force  constants 

on  such  driving-force  dependencies  and  the  degree  to  which  the  experimental 

data  can  be  accounted  for  on  this  basis.  In  addition,  the  electrochemical 

3+/2+  3+/2+ 

activation  parameters  for  Cr(OH_),  and  Eu(OH_)  have  been  evaluated 

zo  z  n 


over  a  wide  range  of  anodic  and  cathodic  overpotentials  and  are  compared  with 
the  corresponding  theoretical  predictions.  Besides  providing  a  more  complete 
description  of  the  driving-force  dependent  reaction  energetics,  such  comparisons 
are  of  particular  interest  since  the  large  (ca.  200  J.  deg  ^  mol 
differences  in  entropy  between  the  oxidized  and  reduced  aquo  cations ^ 
yield  very  substantial  differences  in  the  redox  thermodynamics  at  anodic  and 
ca  tr.ocic  overpotentials . ^  Finally,  some  commonalities  are  explored  between 
the  driving-force  dependence  of  redox  reactivities  for  electrochemical  and 


er.eous  processes , 


t".- 


v  -  •.  *.  •- 


Driving-Force  Dependence  of  Electrochemical  Rate  Constants 

Figures  1  and  2  contain  experimental  Tafel  plots  for  Cr  (OK^)  ^+/^+  and 
3+/?+ 

Eu(OK0)n  ,  respectively,  over  a  range  of  anodic  and  cathodic  overpotentials 

at  the  mercury-aqueous  interface.  These  plots,  which  are  taken  from  ref.  1, 

3+/2+ 

are  shown  as  solid  curves.  Similar  data  have  also  been  obtained  for  VCOH^)^ 

although  over  a  somewhat  smaller  overpotential  range. ^  The  dotted  lines 

in  Figs.  1  and  2  are  the  Tafel  plots  that  would  be  obtained  for  the  same 

value  of  k  at  the  formal  potential,  E-  (i.e.  the  "standard"  rate 
corr  r  f 

constant  kS  )  if  the  transfer  coefficient,  a  ,  equals  0.5  throughout  the 
corr  corr 

entire  overpotential  range.  The  transfer  coefficient  is  defined  as 


a  =  I (RT/F) (dink  /dE)  (1) 

corr  '  corr 

where  the  plus-minus  sign  refers  to  anodic  and  cathodic  rate  constants, 

respectively.  The  strongly  asymmetric  nature  of  the  experimental  anodic  and 

cathodic  Tafel  plots  is  clearly  seen  by  comparison  with  the  dotted  lines. 

The  conventional  expression  for  electrochemical  rate  constants  as  a 

4 

function  of  the  overpotential  can  be  written  for  one-electron  reactions  as 

k  =  A  exp{[-X±F(E-E,)]2MART)  (2) 

corr  i 


where  the  plus-minus  sign  again  refers  to  anodic  and  cathodic  rate  constants, 
A  is  a  preexponential  factor,  and  X  is  the  intrinsic  reorganization  energy. 
This  last  term  is  the  free  energy  required  to  adjust  the  reactant  nuclear 
coordinates  so  that  they  correspond  to  those  of  the  product,  but  without 
electron  transfer.  The  form  of  Eq.  (2)  clearly  predicts  symmetrical  curved 
Tafel  plots  at  anodic  and  cathodic  overpotentials.  This  symmetry  follows 
from  the  assumption  contained  in  Eq .  (2)  that  the  force  constants  of  the 
oxidized  and  reduced  forms  are  equal;  i.e.  the  reactant  and  product 
parabolas  have  identical  shapes. 


. . . 


_  * r  'j.  *_■  : 


We  and  others  have  previously  employed  Eq.  (2)  (or  equivalent  relationships) 


for  comparing  the  rate  constant  -  overpotential  dependence  with  the 

13  12 

harmonic  oscillator  model.  ’  ’  It  is  conventional  to  determine  X  from  the 

experimental  standard  rate  constant  using  Eq.  (2)  and  to  reinsert  this  quantity 

into  Eq .  (2)  to  find  the  predicted  overpotential  dependence  of  log  ^corr 

1  3  12 

which  is  compared  with  the  experimental  data.  ’  ’  A  somewhat  different 

approach  is  used  here  whereby  the  predicted  log  ^corr  “  overpotential 

dependences  are  obtained  directly  from  calculated  free-energy  barriers  derived 

from  the  known  redox  structural  parameters.  This  procedure  is  now  outlined. 

13 

The  work-corrected  rate  constant  can  be  expressed  as 


k  =  Kvr  K  exp (-AG*/RT) 

corr  o  n  n  el 


where  Ko  is  the  statistical  part  of  the  equilibrium  constant  for  forming 

the  precursor  state  from  the  bulk  reactant,  is  the  nuclear  frequency  factor, 

r  is  the  nuclear  tunneling  factor,  <  ..  is  the  electronic  transmission 
n  ex 

coefficient,  and  AG*  is  the  free  energy  of  activation  for  the  elementary 
electron- transfer  step.  Although  the  preexponential  factors  exert 
important  influences  upon  the  absolute  magnitudes  of  ^corr»  the  dependence  of 
kcorr  uPon  the  driving  force  is  expected  to  arise  almost  entirely  from 
variations  in  AG*.  For  the  present  purpose,  it  is  particularly  convenient 
to  formulate  AG*  in  terms  of  intrinsic  reorganization  energies  for  the 
forward  and  reverse  directions,  and  X^,  respectively,  for  the  generalized 
reaction 


Ox  +  e 


Provided  that  the  free  energy-reaction  coordinate  profiles  are  quadratic 
in  nature,  values  of  AG*  as  a  function  of  the  free-energy  driving  force 
_Ge;=  r?(E-E,)]  can  be  obtained  from  the  solution  to  the  following  relations: 


5 


where 


AG*  «  X{X2 

XfX2  =  Xr(l-X)2  +  AGe 


(5a) 

(5b) 


and  X  is  a  dimensionless  parameter  characterizing  the  nuclear  coordinates 
of  the  transition  state^  (see  Fig.  3).  [For  convenience,  quantities  such 
as  X  ,  AG*,  AGe,  etc.  refer  here  and  below  to  the  forward  reaction,  either 
cathodic  or  anodic,  being  considered  even  though  the  reduction  process 
is  commonly  regarded  as  the  "forward  reaction",  as  in  Eq.  (4).]  The 
combination  of  Eqs.  (3)  and  (5)  yields  the  simplified  symmetrical  relation 
Eq.  (2)  only  for  the  special  case  when  Xf  =  Xf  =  X.  However,  when  Xf  J  X^ 
they  instead  predict  nonsymmeirioal  log  ^corr  -  overpotential  plots  in  the 
anodic  and  cathodic  directions.  The  latter  will  clearly  apply 
when  the  force  constants  differ  in  the  oxidized  and  reduced  forms. 

The  various  components  of  Eqs.  (3)  and  (5)  can  be  evaluated  (or  estimated) 
as  follows.  The  forward  and  reverse  reorganization  energies  for  Eq.  (4) 
can  be  separated  into  inner -shell  (reactant  bond  distortional)  and  outer- 
shell  (solvent  reorganization)  contributions  according  to 


X 


f 


X 

os 


X 

r 


X 

os 


(6a) 

(6b) 


where  X  and  X^g  are  the  inner-shell  reorganization  energies  associated  with 
the  reactant  and  product  species,  respectively. 

These  may  be  related  to  the  corresponding  force  constant ,'  f^,  for 
symmetrical  bond  stretching  (or  compression)  of  the  ith  bond,  and'  the  change 
in  bond  distance  between  the  oxidized  and  reduced  forms,  by  means  of 


where  v  is  the  bond  vibrational  frequency,  c  is  the  velocity  of  light,  and 
is  the  reduced  mass  of  the  vibrator. 

The  outer-shell  reorganization  energy,  A  ,  will  be  assumed  to  be  equal 
in  the  oxidized  and  reduced  forms  {.vide  infra).  It  is  conventionally  given  by 


where  e  is  the  electronic  charge,  N  is  Avogadro's  number,  r  is  the  reactant 

radius,  R  is  the  distance  from  the  reactant  to  its  image  in  the  metal 

surface,  and  e  and  c  are  the  optical  and  static  dielectric  constants, 
op  s 

o  in 

respectively.  For  the  present  aquo  reactants,  r  was  taken  as  3.25  A 
„  o  16 

and  R  as  13  A. 

3+  /  2«f 

Values  of  AG*  for  the  CrCOE^)^  couple  were  calculated  as  a  function 

of  overpotential  using  Eqs.  (5)  -  (8),  with  p  taken  as  the  mass  of  the  aquo 


ligand,  along  with  the  experimental  values  vm  =  543  cm  vJ1  =  380  cm 
and  la  ■  0.20  This  yields  inner-shell  reorganization  energies,  A^  and 

T  *1 

equal  to  223  and  110  kJ  mol  ,  respectively;  and  overall  reorganization 
energies,  A^  and  A^_,  of  310  and  197  kJ  mol  ,  respectively.  The  corresponding 
calculated  values  of  kcorr  were  obtained  by  combining  the  values  of  AG* 
with  numerical  estimates  of  the  preexponential  factors  in  Eq.  (3).  Details 
concerning  the  physical  significance  of  these  quantities  are 


-1  9,17 


given 


13  18  19  o 

elsewhere.  ’  *  The  composite  quantity  was  taken  as  0.5  A; 


similar  values  (within  ca.  twofold)  are  indicated  from  both  experimental  and 


theoretical  lines  of  evidence.  The  nuclear  tunneling  factor  T  was  obtained 

21 

from  the  driving-force  dependent  relationship  of  Holstein  as  expressed 

22 

conveniently  by  Sutin.  It  is  calculated  to  equal  3.4  at  E=E^,  and  depends 

only  slightly  upon  overpotential  in  the  range  of  interest  here.  The  nuclear 

13  -1 

frecuencv  factor  .  was  determined  as  1.0  x  10  sec  from  the  above  stretch! 
n 


encies  as  described  in  ref.  13. 


The  resulting  calculated  values  of  log  k^orr 


3+/2+ 


against  (E-E^)  for  CrCOH^)^ 
in  both  anodic  and  cathodic  directions  are  plotted  as  dashed  curves  in  Fig.  1. 

‘  5~ 

The  standard  rate  constant,  k  ,  obtained  from  theorv  is  in  almost  exact 

corr 

agreement  with  the  experimental  value,  3  x  10-^  cm  sec~\  Given  the 

uncertainties  and  assumptions  involved  in  deriving  the  former  (probably 

ca.  ±20-fold)  this  agreement  is  somewhat  fortuitous.  For  convenience, 

the  calculated  and  experimental  values  of  k~orr  are  set  exactly  equal  in 

Fig.  1  so  to  facilitate  the  comparison  of  the  theoretical  predictions  for 

the  driving-force  dependence  of  kcorr  with  experiment. 

The  theoretical  analysis  for  Eu (0^)  cannot  be  performed  on  the 

same  quantitative  basis  as  for  Cr (0H2) ^+/2+  since  experimental  data  for 

v  ,  v__,  and  Aa  are  unavailable.  However,  approximate  calculations  of  the 

log  “  overpotential  plots  for  this  couple  may  still  be  carried  out. 

s  ^5  “1 

From  the  experimental  standard  rate  constant,  k  =  3  x  10  cm  sec  ,  together 

corr 

with  the  theoretical  value  of  A^  assuming  r  *  3.5  i,  R  =  13  %  [87  kj.  mol~\ 

Eq.  (8)],  and  noting  that  the  reaction  is  likely  to  be  nearly  adiabatic 
23 

'1),  an  intrinsic  inner-shell  barrier,  AG*nt  (i.e.  AG*  for  AGe  =  0) 

of  ca.  30  kj.  mol  ^  is  inferred.  Since  the  interaction  of  Eu(III)  and 

Eu(II)  with  the  aquo  ligands  is  predominantly  electrostatic,  the  Eu-OHj 

force  constants  should  be  approximately  proportional  to  the  ionic  charge 

24 

divided  by  the  cation-aquo  bond  distance,  a.  Given  that  a  is  likely 
o 

to  be  about  2  A,  and  the  difference  in  bond  lengths  between  Eu(III)-0H2  and 

o  23 

Eu(lI)-0Ko  is  around  0.22  A,  from  the  above  values  of  AG*  and  A  and 
2  mt  os 

utilizing  Eq.  (5)  for  AGe  =  0  along  with  Eq.  (6)  leads  to  estimates  of 

Ox  Red  -1 

A."  and  A,  of  170  and  95  kj.  mol  ,  respectivelv .  The  overpotential 
in  in  t  r  . 

dependence  of  log  k  for  Eu(0H„)  ^+/^+  uredicted  bv  these  parameters 
r  s  corr  2  n  -  -  r 

are  shown  as  dashed  curves  in  Fig.  2. 


Comparison  between  the  experimental  and  calculated  Tafel  plots  for 

Cr (OH.)  and  Eu(OH0)  •3+//2+  in  Figs.  1  and  2  reveals  that  the  observed 

2  o  2  n 

anocic-cathodic  asymmetry  is  accounted  for  in  part  by  the  inclusion  of 
larger  force  constants  for  the  oxidized  compared  to  reduced  forms  in  the 
theoretical  calculations.  Thus  these  ratios  of  inner-shell  force  constants, 
estimated  to  be  2.05  and  ca.  1.8  for  Cr (0H„)  and  Eu(OH0)  2+/,2+, 

respectively,  yield  noticeably  smaller  slopes  in  the  Tafel  plots  at 
anodic  as  compared  to  cathodic  overpotentials.  Even  though  the  predicted 
Tafel  plots  at  cathodic  overpotentials  do  exhibit  some  curvature, 
these  are  nevertheless  closely  similar  to  the  essentially  linear 
experimental  plots  in  Figs.  1  and  2.  At  anodic  overpotentials,  the  marked 


curvature  in  the  experimental  Tafel  plots  is  roughly  mimiced  by  the 

2+ 

calculated  plots,  especially  for  C^OH^)^  oxidation. 

A  comparison  between  the  experimental  and  calculated  transfer  coefficients 


for  Cr  (OH  .,)  ,3+/,2+,  derived  from  the  slopes  of  the  curves  in  Fig.  1  [Eq.  (1)], 
Z  0 

£ 

is  given  in  Table  I.  The  calculated  anodic  and  cathodic  quantities,  ccalc 
and  aca2c  respectively,  were  obtained  from  the  dashed  curves,  i.e.,  using 


Q.  C 

individual  force  constants.  Although  the  changes  in  &ca-]_c  and  acaic  wit*1 


varying  overpotential  are  similar,  the  former  are  substantially  larger  at 


a  given  overpotential,  in  harmony  with  the  behavior  of  the  corresponding 

experimental  quantities  °ca^c  an^  Qcaic  *) •  In  contrast,  the  transfer 

coefficient,  aS^f  ,  calculated  from  the  conventional  treatment  [Eq.  (2)] 
calc 

displays  a  symmetrical  overpotential  dependence,  in  marked  disagreement 

X  3 

with  the  observed  behavior  (Table  I) .  The  values  of  a  are  nevertheless 

corr 

g 

significantly  smaller  than  a  1  at  large  anodic  overpotentials;  this  is 


also  evident  from  the  deivations  between  the  solid  and  dashed  curves  in 


These  residual  discrepancies  betveen  theory  and  experiment  at  anodic 

over potentials  may  be  due  to  several  causes.  Most  simply,  the  actual 

curvature  in  the  experimental  log  k  -E  plots  nav  be  somewhat  different  than 

corr 

that  shown  in  Figs.  1  and  2  due  to  possible  potential-dependent  errors  in  the 
application  of  double-layer  corrections  to  the  observed  rate  constants. 
Although  the  extent  of  this  uncertainty  is  not  large  (around  ca.  twofold 
in  log  ^corr) »  is  greatest  at  the  largest  anodic  overpotentials  in  view 
of  the  substantial  anionic  specific  adsorption  encountered  even  in  hexafluoro- 
phosphate  electrolytes.1  However,  it  is  more  likely  that  the  discrepancies 
are  due  largely  to  remaining  inadequacies  of  the  theoretical  model.  One 
possibility  is  that  the  metal-ligand  vibrational  modes  are  anharmonic. 

However,  trial  calculations  using  anharmonic  Morse  potentials  in  place  of  the 


harmonic  oscillator  model  [Eq.  (7)]  yielded  almost  imperceptible  changes  in 
tne  Tafel  plots,  at  least  for  the  conditions  encountered 

2  5  o  «  /  24- 

in  Fig.  1.  A  more  important  factor,  at  least  for  Cr(OH-),  '  'r  is  liable 

Z  b 

to  reside  in  the  above  assumption  that  all  the  metal-ligand  bonds  are 

2+ 

equivalent.  A  pair  of  tvzns  chromium-oxygen  bonds  in  CrCOH^)^  are  expected 
to  be  weaker,  and  therefore  have  smaller  force  constants,  than  the  other 
four  due  to  Jahn-Teller  distortion.  Since  these  bonds  are  also  especially 
elongated, ^  they  may  well  contribute  importantly  to  the  inner-shell 

barrier,  yielding  a  more  pronounced  curvature  in  the  calculated  anodic 
Tafel  plot. 

As  noted  above,  we  previously  attributed  the  observed  Tafel  plot 

1  3 

asymmetry  chiefly  to  the  influence  of  short-range  solvent  polarization.  * 

Although  a  large  fraction  of  the  asymmetry  effect  now  appears  to  reside  in 

the  inner-shell  barrier,  the  present  analysis  suggests  a  way  in  which  the 

outer-shell  reorganization  can  contribute  as  well.  We  assumed  above 

that  a  is  equal  for  the  reduction  and  oxidation  directions:  i.e.  the  outer - 
os 

shell  "force  constant"  is  independent  of  the  metal  oxidation  state.  This 

follows  from  the  simple  dielectric  continuum  model  used  to  estimate  the 

outer-shell  barrier.  If  reorganization  of  hydrogen -bonded  outer -shell 

water  molecules  contributes  signficantly  to  the  free-energy  barrier,  greater 

Tafel  plot  asymmetry  should  result  since  the  effective  force  constants 

associated  with  such  motions  should  also  be  larger  in  the  higher  oxidation 
3 

state.  However,  quantitative  calculations  of  such  effects  are  precluded 
at  present  due  to  the  incompleteness  of  molecular  structural  information 
for  the  surrounding  solvent. 

26  27 

With  regard  to  recent  provocative  discussions,  ’  we  stress  that 
it  is  necessary  to  include  the  outer-shell  solvent  contribution  to  the 


reorganization  barrier  in  order  to  obtain  satisfactory  agreement  between 

the  theoretical  and  experimental  rate  constants  for  Cr  (OH,,) , 3+^2*  This  is  also 

2  6 

the  case  tor  a  variety  of  other  homogeneous  and  electrochemical  reactions.1^ 

To  illustrate  this  point  for  Cr (0H-)  3+/,2+,  Fig.  1  also  contains  log  k  -E 
plots  calculated  by  ignoring  the  solvent  contribution  (dotted-dashed  lines). 
Although  exhibiting  very  similar  dependences  of  log  kcorr  upon  overpotential 
to  that  seen  with  the  proper  inclusion  of  the  solvent  reorganization  component,  the 
calculated  values  of  log  kcorr  are  substantially  (ca.  10A  fold)  larger  than 
the  experimental  values. 

In  view  of  the  continuing  interest  and  speculation  concerning  the 

virtues  of  Tafel  plots  for  examining  some  fundamental  aspects  of  electro- 
26  28  29 

chemical  kinetics,  1  it  is  worthwhile  to  note  some  more  general  implications 

of  the  foregoing.  It  is  evident  that  the  degree  of  curvature  of  the  Tafel 

plots  can  be  at  least  as  sensitive  to  the  degree  of  mismatch  between  the 

forward  and  reverse  reorganization  energies  as  to  the  absolute  magnitude 

of  these  energies,  and  hence  to  the  intrinsic  barrier,  AG*  ,  and  kS 

int’  corr 

It  is  therefore  extremely  dangerous  to  estimate  intrinsic  barriers  from  the 

28,30,31  „  t 

curvature  of  Tafel  plots.  Further,  the  Tafel  plot  shapes  are  not  very 

diagnostic  of  the  factors  contributing  to  reactant  activation,  such  as  the 

role  of  solvent  reorganization.  A  more  valuable  approach  to  examining  such 

fundamental  questions  is  to  compare  the  absolute  magnitudes  of  the  calculated 

19 

and  experimental  rate  constants  under  well-defined  conditions  rather  than 
their  dependence  upon  the  thermodynamic  driving  force.  Nevertheless,  such 
Tafel  comparisons  do  provide  one  means  of  evaluating  more  subtle  driving-force 
dependent  effects  with  the  aim  of  "fine  tuning"  theoretical  formulations. 


11 


Activation  parameters  clearly  provide  an  additional  measure  of  information 

on  the  reaction  energetics  beyond  that  yielded  by  rate  measurements  at  a 

single  temperature.  These  parameters  can  be  particularly  revealing,  for 

example,  of  the  solvational  changes  accompanying  'transition-state  formation 

32 

for  simple  redox  reactions.  Two  distinct  types  of  activation  parameters 

can  usefully  be  distinguished  for  electrochemical  reactions.  The  first 

type,  the  so-called  "real"  activation  parameters,  refer  to  the  temperature 

33 

dependence  of  the  rate  constant  measured  at  a  fixed  overpotential . 

The  second  type,  so-called  "ideal"  activation  parameters,  are  obtained  from 

the  temperature  derivative  evaluated  at  a  fixed  nonisothermal  electrode 
32-34 

potential.  The  difference  between  these  two  types  lies  in  the  thermo¬ 

dynamic  entropy  change,  the  "reaction  entropy"  AS°c>  brought  about  by  electron 
transfer,  i.e.  the  entropy  difference  between  the  reduced  and  oxidized 
forms  of  the  redox  couple.'1'1  This  entropy  change  is  associated  with  the 
temperature  dependence  of  the  standard  electrode  potential,11  and  gives 
rise  to  an  inherent  asymmetry  of  the  electrochemical  reaction  energetics 
in  the  anodic  and  cathodic  directions.  Thus  at  the  standard  potential 
where  the  free-energy  driving  force,  AG°,  necessarily  equals  zero,  there 
will  nevertheless  be  a  nonzero  (and  often  substantial)  reductive  enthalpic 

driving  force,  AK°  ,  since  then  AH°  *  TAS°  .^2  ^  Thus  for  Cr(OH_) 

°  ’  re’  rc  rc  2  6  * 

AS°  =  205  J.  K  1  mol  1,11  and  therefore  at  Er,  AH°  *  61  kj.  mol'1.  This 
rc  f  rc 

is  a  key  feature  which  sets  apart  electrochemical  exchange  reactions  from 

the  corresponding  homogeneous  self-exchange  reactions;  the  latter  are 

inherently  symmetrical  since  they  can  be  conceived  as  a  coupled  pair  of 

2 

electrocnemical  exchange  reactions. 


.■  ././v  a.’ 


The  "ideal"  activation  parameters  can  be  identified  with  the  actual 
enthalpic  and  entropic  barriers  at  the  electrode  potential  at  which  they 
are  evaluated.  The  cathodic  and  anodic  "ideal"  parameters  will  therefore 
differ  markedly  even  at  the  standard  electrode  potential.  The  "real" 
activation  parameters,  on  the  other  hand,  contain  a  correction  for  the 
entropic  driving  force;  they  will  necessarily  be  identical  in  the  anodic 
and  cathodic  directions  at  the  standard  potential  since  the  temperature- 
dependent  anodic  and  cathodic  rate  constants  will  be  equal  under  these 
conditions.  The  latter  are  therefore  closely  related  to  the  activation 
parameters  for  homogeneous  self -exchange  reactions.  However,  the  "ideal" 
parameters  are  more  appropriate  for  the  present  purposes  since  they  reveal 
the  true  asymmetric  nature  of  the  electrochemical  reaction  energetics. 

Table  II contains  a  comparison  of  experimental  "ideal"  enthalpies  and 

entropies  of  activation,  AH*  and  AS*  ,  respectively,  with  the 
c  corr  corr 

corresponding  calculated  quantities,  AH*a^c  and  dS*alc,  The  work-corrected 
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activation  enthalpies,  AH*  ,  were  obtained  from 

corr 

“Jorr  '  -*B1"kcorr/8(1/T>¥i 

where  the  temperature  dependence  of  k  is  evaluated  at  the  constant 

corr 

nonisothermal  cell  potential  En^.  The  values  of  k  were  evaluated  over 

r  corr 

the  temperature  range  ca.  5—50 0 C ;  experimental  details  are  given  in  ref.  32. 

At  cathodic  overpotentials  the  kcorr  values  were  obtained  from  apparent 

rate  constants  measured  in  1  M  NaClO^  (pH  2.5)  as  outlined  in  ref.  32. 

At  anodic  overpotentials,  they  were  evaluated  from  rate  data  measured  in 

0.4  M  KPF,  (pH  2.5)  as  outlined  in  ref.  1.  [Since  quantitative  double-laye 

—  D 

data  are  not  available  for  0.4  K  KPF^  at  potentials  positive  of  the  point 

of  zero  charge  (where  Pr&”  specific  adsorption  occurs),  the 

r'ou'-'e- laver  corrections  are  assumed  to  be  independent  of  temperature.  The 


approximate  validity  of  this  assumption  is  supported  by  a  similar  finding  for 
3+  32 

Cr ( OK^ ) £  reduction  “  and  the  small  absolute  magnitude  of  the 

1 

corrections  required  for  Cr (0K2) ^  oxidation  in  0.4  M  KPFfe.  ]  The  work- 

corrected  entropies  of  activation,  AS£o„r,  also  listed  in  Table  II,  were 

obtained  from  the  corresponding  values  of  AH*  and  k  at  a  given 

corr  corr  6 

temperature,  along  with  the  above  estimates  of  K  <  ,  and  v  using 

o  el  n  6 


AS*  =  Rink  -  Rln(K  k  v  )  +  AH*  /T 
corr  corr  o  el  n  corr 


(10) 


Kote  that  the  temperature-dependent  nuclear-tunneling  factor  F  is  contained 

n 

within  -H*orr  and  AS£orr>  so  that  the  latter  can  be  regarded  as  "semiclassical" 

.  .  22 
quantities . 

The  corresponding  calculated  quantities,  AH*  .  and  AS*  ,  ,  were  obtained 

calc  calc 

follows.  As  shown  in  the  Appendix,  differentiation  of  Eq.  (5)  with  respect 
to  temperature  yields  an  expression  for  the  classical  activation  entropy 
(i.e.  not  containing  the  nuclear-tunneling  correction): 

AS*  =  AS' A^X  [  A2  +  (Af  -  Ar)(AG°  +  A^.)]"172  (11) 
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by  assuming  that  dAf/dT  *  dAr/dT  •=  0,  where  AS'  is  the  entropic  driving 
force  in  the  cathodic  or  anodic  direction,  as  appropriate.  An  equivalent 
expression  for  AS*  is  simply  (see  Appendix) 

AS*  =  aAS°  (12a) 

where  the  transfer  coefficient,  a,  as  operationally  defined  by  Eq.  (1)  is 
calculated  from 

a  =  AfX(XAf  +  Ar  -  XAr)_1  (12b) 

Although  Eqs.  (11)  and  (12a)  are  equivalent,  the  latter  formulation  is  more 
revealing  since  it  separates  the  simple  phenomenological  definition  of  IS* 


[Eq.  (12a)]  and  its  detailed  calculation  based  on  a  harmonic  oscillator 
treatment  for  a  [Eq.  (12b)]. 


The  classical  activation  enthalpy,  AH* ,  is  obtained  from  AS*  together 
with  the  classical  activation  free  energy  from  Eq.  (5)  using 

AH*  =  AG*  +  TAS*  (13) 

The  calculated  value  of  AH*.  AH£alc>  desired  for  comparison 

with  the  experimental  value,  AH*orr, is  obtained  by  including  the  temperature- 
dependent  nuclear  tunneling  correction: 

AH*  .  *  AH*  -  R[dlnr  /d(l/T)]  (14) 

calc  n 

Similarly,  the  corresponding  calculated  activation  entropy  can  be  obtained 
from 


AS*  .  -  AS*  -  Rlnl  -  (R/T)dlnT  /d(l/T)  (15) 

calc  n  n 


The  resulting  values  of  ^H*a^c  and  for  CrCOI^)^*^*  are  listed 

alongside  the  corresponding  experimental  estimates  in  Table  II*  The  nuclear  - 


tunneling  corrections  to  AH£alc  and  AS*alc  range  from  ce.  -3.5  to  -10 
kj.  mol  ^  and  -7  to  -22  J.  K  ^  mol  \  respectively,  becoming  more  negative 
with  increasing  .negative  .electrode  potential. 


Inspection  of  Table  II  reveals  that  the  striking  numerical  differences  in 

the  calculated  activation  parameters  for  the  anodic  and  catholic  directions 

brought  about  by  the  large  reaction  entropy  for  (^(OH.^)^  is  matched 

to  a  large  degree  by  the  experimental  values.  This  agreement  extends,  albeit 

2+ 

approximately,  to  the  very  small  values  of  AH*  and  AH*  ,  for  Cr(0Ho), 
rr  '  ’  J  corr  calc  2  6 

oxidation  at  large  anodic  overpotentials  and  their  mild  overpotential 

3+ 

dependence  in  comparison  with  that  for  CrtO^)^  reduction  at  cathodic 


■  -  'W'ifw) 


15 


overpotentials.  Similar  calculations  for  Eu(OH-) 


3+/ 2+ 


are  not  feasible 


rpOlciU.XCi.9  .  uoj.  w  t  x  X  w  *.  i-U*  \  ^ 

since  the  required  quantitative  structural  data  unavailable  for  this  couple. 

3+/ 2+ 

Nevertheless,  the  experimental  values  of  IK*  and  AS*  for  Eu(OH-) 

’  r  corr  corr  2  n 

also  display  striking  differences  in  the  anodic  and  cathodic  directions  as 

a  consequence  of  the  large  AS°c  value  (195  J.  K  ^  mol  ^  for  this  couple. 

In  addition,  similarly  small  and  even  negative  values  of  AH*  for 

corr 

2+ 

Eu(OH-)  oxidation  were  obtained  at  anodic  overpotentials  [eg.  AH*  = 

2  n  corr 

-5  kJ.  mol-1  at  (E-Ef )  =  400  mV]. 

An  apparently  unexpected  finding  which  we  alluded  to  earlier^"  is  the 

observation  of  near  zero  activation  enthalpies  at  anodic  overpotentials 

even  though  the  activation  entropies  remain  large  and  negative  under  these 

conditions  (Table  II).  Intuitively,  one  might  expect  the  activation  entropy 

to  equal  zero  in  the  absence  of  an  enthalpic  barrier.^  In  order  to  examine 

this  matter  further,  values  of  ^G*a^c»  ^*aic*  an(*  ^*alc  were  calculated  for 
2+ 

Cr(0Eo),  oxidation  over  a  considerable  range  of  free  energy  driving 

<£  O 

forces,  AG°,  corresponding  to  a  much  greater  span  of  anodic  overpotentials 

than  can  be  examined  experimentally.  The  resulting  plots  of  ^G*a^c»  ^*aic> 

and  4S*£^c  against  AG®  are  given  in  Fig.  4.  (Nuclear  tunneling  corrections 

were  omitted  since  these  are  small  and  the  approximations  inherent  in  their 

38 

calculation  become  somewhat  ambiguous  for  highly  exoergonic  reactions.  ) 

Figure  4  indicates  that  £H*a2c  decreases  sharply  and  even  assumes  small  negative 

values  at  substantially  smaller  driving  forces  than  are  necessary  for  4G*a^c 

to  become  zero.  Indeed,  AH*a^c  attains  a  value  of  zero  at  two  points, 

approximately  when  AG°  *=  -(X^  +  2TAS°)  and  when  AGe  *  ~/'ft  whereas  4G*_jc 

and  AS*  .  become  zero  onlv  when  AG°  =  -X, .  The  significance  of  negative 
calc  '  I 

activation  enthalpies  for  related  homogeneous  reactions  has  been  discussed 


Marcus  and  Sutin. 
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A  related  feature  of  Fig.  4  is  the  appearance  of  an  enthalpic  "inverted 

region",  i.e.  a  realm  of  increasing  activation  enthalpy  with  increasing 

exoergonicity  within  the  "normal  free-energv  region",  i.e.  where  the  activation 

free  energy  decreases  with  increasing  driving  force.  It  has  been 

40 

noted  that  the  much-discussed  free-energy  inverted  region  should  be  absent  for 

41  42 

electrochemical  reactions  at  metal  surfaces.  ’  This  is  because  oxidative  . 
electron  transfer  to  metal  energy  states  above  the  Fermi  level,  or  reductive 
electron  transfer  from  states  below  the  Fermi  level,  will  provide  an 
activationless  pathway  even  when  -AG®  exceeds  A^ .  However,  the  enthalpic 
inverted  region  should  nevertheless  be  observable  for  electrochemical 
reactions.  It  would  be  interesting  to  explore  this  possibility  experimentally  by 
selecting  reactions  with  sufficiently  negative  values  of  AS®  so  that  the 
anticipated  enthalpic  inverted  region  occurs  at  sufficiently  large  values 
of  AG*orj.  corresponding  to  measurably  slow  reaction  rates. 

It  should  be  noted  that  the  anodic-cathodic  driving  force  asymmetry  seen 
for  the  activation  parameters  has  a  different  origin  than  that  of  the  rate 
constants  discussed  above.  Thus  purely  symmetrical  Tafel  plots  will  be 
predicted  provided  that  A^  =  A^  irrespective  of  major  differences  in  the 
entropic  and  enthalpic  components  of  these  reorganization  energies.^-  Nevertheless, 
these  two  types  of  driving  force  asymmetry  are  in  a  sense  related  since  they 
both  arise  from  equilibrium  structural  changes  brought  about  by  electron 
transfer.  In  addition,  we  have  speculated  that  the  especially  small  transfer 
coefficients  observed  at  large  anodic  overpotentials  (Figs.  1,2)  may  be 
connected  with  the  occurrence  of  an  enthalpically  barrierless  process 
(i.e.  AH*oi^  -  0)  under  these  conditions.^  Although  the  above  analysis  based 
on  force  constant  differences  largely  removes  the  need  to  invoke  such 
"anomolous  activation"  mechanisms,  such  notions  should  nonetheless  not  be 


wholly  discarded  at  this  point. 


The  foregoing  demonstrates  that  the  chief  property  of  simple  electro¬ 
chemical  reactions  which  leads  to  unusual  driving  force-dependent  reactivities 
is  their  inherent  chemical  asymmetry.  In  this  regard  interfacial  redox  reactions 
have  much  in  common  with  homogeneous  cross  reactions  between  structurally 

diverse  redox  centers,  and  in  fact  are  more  akin  to  these  than  to  homogeneous 

3+/2+ 

self-exchange  reactions.  As  an  illustration,  consider  the  V(OH2>6 

3+/2+ 

couple.  The  VCOH^)^  homogeneous  self-exchange  reaction  is  a  completely 

symmetrical  process  and  is  therefore  necessarily  isoenergetic  as  well  as  thermo¬ 
neutral  and  isoentropic.  In  contrast,  V(OH2)^+  electrochemical  exchange  is 
characterized  by  net  AS®  and  AH®  values  of  -155  J  K  1  mol  1  and  -46  kj  mol  1, 
respectively,11  for  the  oxidation  reaction  at  the  standard  potential  (i.e.  for 
1G°  =  0) .  These  thermodynamic  changes  are  closely  shared  by  the 
oxidation  of  V(0H2)62+  by  Cr(bpy)33+  in  homogeneous  aqueous  solution  (bpy  = 

2,2 ’-bipyridine) ,  since  the  formal  potentials  of  the  reacting  redox  couples 
are  virtually  identical  (-480  mV  vs.  s.c.e.11’^3)  ,  and  AS®c  for  Cr(bpy)3  ^ 
is  close  to  zero  (15  J.  K  mol"  .  [The  values  of  AS®  and  AH®  are  deduced 
to  be  -140  J.  K-1  mol-1  and  -41  kj.  mol”1,  respectively'.]  It  is  not  surprising, 
therefore,  that  cross  reactions  between  a^o  and  low-spin  polypyridine  complexes 
exhibit  similarly  "unusual"  activation  parameters  to  electrochemical 
reactions  involving  the  former  reactants.3  For  example,  the  highly  exoergonic 
oxidations  of  Fe(OH2)g)2+  by  Ru(bpy)33+,  Ru(phen)3  and  Fe(phen)3  are 
characterized  by  small  negative  activation  enthalpies  and  large  negative 

activation  entropies. ^  3  It  has  been  shown  that  the  substantial  negative 

-1  -1  12 

entropic  driving  force  for  these  reactions  (AS®  -  -175  J  K  mol  )  can 

account  in  broad  terms  for  this  kinetic  behavior  on  the  basis  of  the  measured 

39 

rate  parameters  for  the  parent  self-exchange  reactions. 


.  L'  LfrY 


18 


Moreover,  comparably  "anomolous"  rate-driving  force  dependencies  to  those 

discussed  here  for  electrochemical  oxidations  have  also  been  observed  for  a 

2 

number  of  homogeneous  oxidations  of  aquo  complexes.  Analysis  of  the  driving- 

force  dependencies  of  homogeneous  reactivity  is  complicated  by  the  need 

to  account  for  activation  of  the  reducing,  as  well  as  oxidizing,  redox  center. 

Low-spin  polypyridine  complexes  provide  particularly  useful  coreactants  in 
2  43 

this  regard  since  the  small  entropic  changes  and  intrinsic  activation 
47 

barriers  that  are  characteristic  of  such  couples  enable  the  reaction 

energetics  to  be  dominated  by  the  other  redox  center.  In  this  respect 

such  polypyridine  reactants  provide  similar  redox  characteristics  to  metal 

surfaces.  Although  the  "overpotential"  provided  by  the  former  is  fixed  by 

its  formal  potential,  driving  force  dependencies  can  nevertheless  be  scrutinized 

by  selecting  a  related  series  of  polypyridine  oxidants  having  a  suitable 

44 

range  of  formal  potentials. 

The  dependence  of  homogeneous  reactivities  upon  the  driving  force  has 
2  5  48  49 

often  been  discussed  ’  ’  ’  in  terms  of  observed  deviations  from  the 

4 

"Marcus  cross  relationship" 


^12  '  (kllk22K12f) 


(16a) 


where 


f  =  (log  K12)2/[4  log(k11k22/A2) ]  (16b) 


Here  k^2  is  the  rate  constant  for  a  cross  reaction  having  an  equilibrium 
constant  K^2  from  the  known  self -exchange  rate  constants,  k^  and  k22,  for 
the  constituent  redox  couples.  Since  Eq.  (16)  is  derived  by  assuming  that 
the  force  constants  for  the  oxidized  and  reduced  form  of  both  reacting 
couples  are  equal,  ’  progressively  larger  deviations  from  this  relation 
are  expected  as  the  driving  force  is  increased.  When  (as  for  aquo 


complex  oxidations) ,  the  observed  values  of  £re  expected  to  be  progressively 
smaller  than  those  predicted  from  Eq.  (16)  as  increases. ^  Such  considerations 
are  entirely  analogous  to  those  considered  above  for  electrochemical  reactions. 

Thus  the  symmetrical  Tafel  relations  obtained  from  the  harmonic  oscillator  model 
with  [Eq.  (12)]  can  be  derived  from  Eq.  (16)  by  noting  that  a  metal 

surface  can  be  regarded  formally  as  a  coreactant  with  variable  driving  force 
and  zero  intrinsic  barrier;  (i.e.  having  a  "self-exchange  rate  constant" 
equal  to  the  preexponential  factor,  A). 

Chou  et  al  regarded  the  deviations  from  Eq.  (16)  due  to  unequal  force 
constants  as  being  too  small  to  account  for  the  observed  discrepancies  for 
a  number  of  homogeneous  cross  reactions,  including  several  oxidations  of 


acuo  complexes.  Undoubtedly  other  factors,  such  as  additional  components 

.....  5  48  49 


of  work  terms,  nonadiabaticity ,' 


use  of  inappropriate  self-exchange 


rate  constants,  etc.  also  contribute  importantly  in  many  cases.  Nevertheless, 
the  above  analysis  clearly  indicates  that  force  constant  differences  can 
provide  significant  driving  force-dependent  deviations  from  Eq.  (2)  and 
therefore  from  Eq.  (16).  Indeed,  by  employing  the  free-energy  form  of  Eq.  (16) 
we  have  demonstrated  that  the  observed  driving  force-dependent  deviations 
observed  for  homogeneous  oxidations  of  aquo  complexes  are  quantitatively 
consistent  with  the  corresponding  deviations  for  the  electrooxidation  reactions 
that  are  discussed  here.^ 

The  present  analysis  focusses  attention  on  aquo  redox  couples  since 
structural  as  well  as  extensive  rate  parameter-overpotential  data  are 
available  for  these  systems.  Nevertheless,  such  rate-driving  force  asymmetries 
should  be  observed  for  any  redox  couple  with  suitably  large  force  constant 
differences  between  the  oxidized  and  reduced  forms,  providing  that  inner-shell 


reorganization  provides  an  important  contribution  to  the  activation  barrier. 


These  two  factors  will  tend  to  occur  together  since  it  is  anticipated  that 
redox  couples  having  such  large  force  constant  differences  will  also  exhibit 
sizable  bond-length  changes  upon  electron  transfer,  and  hence  have  large 
inner-shell  barriers.  Since  such  large  structural  changes  commonly  give 
rise  to  rapid  irreversible  chemical  steps  following  electron  transfer, 
the  number  of  chemically  reversible  couples  having  the  required  properties 
may  be  relatively  small.  Nevertheless,  the  possible  occurrence  of  such  Tafel 
plot  asymmetry,  yielding  otherwise  unexpected  Tafel  linearity  or  curvature, 
should  be  borne  in  mind  when  examining  the  electrochemical  kinetics  of  reactions 
involving  coupled  chemical  steps,  including  those  involving  proton  transfer 
such  as  proton  electroreduction. 
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APPENDIX:  Equivalence  of  Equations  11  and  12 


The  free  energy  expression 

AG*  =  X2  X 

nay  be  differentiated  with  respect  to  temperature  to  yield 

AS*  =  -2Af  X  dX/dT 

By  solving  Eq.  5  one  obtains  for  X: 

x  =  {_xr  +  [X2  +  ( Af  -  Xr)(AGe  +  Xr)]1/2}  (Xf  -  Xr)-1 

Assuming  that  dX/dT  *  0,  differentiation  of  Eq.  A2  yields: 

dX/dT  =  -  \  AS°[X2  +  (X  -  X  ) (AG°  +  X  )]"1/2 
2  r  f  r  r 

Substitution  of  Eq.  A3  into  Eq.  Al  yields  Eq.  11: 

AS*  =  AS°XX.[X2  +  (X  -  X  ) (AG°  +  X  ) ]- 1/2 
f  r  f  r  r 

Rearranging  Eq.  A2  we  find  that: 

[X2  +  (X,  -  X  )(AG°  +  X  )]1/2  =  X(X,  -  X  )  +  X 
r  f  r  r  f  r  r 

and 

[X2  +  (X,  -  X  )  (AGe  +  X  )]_1/2  =  (XX,  +  X  -  XX  r1 
r  f  r  r' J  f  r  r 

Incorporating  this  result  into  Eq.  11  yields: 

AS*  =  AS0  X,  X(XXr  +  X  -  XX  )~1 
f  f  r  r 

From  Eq.  5b  one  can  obtain  (noting  that  dX/dX  =  0): 

diG6/ d>:  =  2  (:•::•  +  x_  -  xxj  (a7) 


From  Eq.  5a  one  finds  that 


dAG*/dX  =  2X  A  (A8) 

In  view  of  Eqs.  A  7  and  A  8,  Eq.  A.6  may  be  rewritten  as 

AS*  =  ASe(dAG*/dX)(dAG0/dX)_1  *  ASe  dAG*/dAG°  (A9) 

Noting  that  the  transfer  coefficient  a  is  defined  as  dAG*/dAG6,  we  conclude 
that 

AS*  =  aAS°  (12a) 

At  first  sight,  Eq.  (12a)  appears  to  be  simply  the  entropic  analog 
of  the  free-energy  expression 

AG*  =  AG*  +  aAG°  (A10) 

int 

where  AG*^  is  the  "intrinsic"  free-energy  barrier,  i.e.  that  which  remains  ir 

the  absence  of  the  free-energy  driving  force  AG°.  [Equation  (A10)  is  the  most 

general  form  of  the  dependence  of  the  activation  barrier  upon  the  driving 

force.]  Indeed,  the  above  assumption,  dX/dT  =  0, made  in  deriving  Eq.  (12a) 

is  equivalent  to  asserting  that  the  intrinsic  entropic  barrier  equals  zero. 

However,  a  difference  between  the  forms  of  Eqs.  (12a)  and  (A10)  is  that  the 

transfer  coefficient  appearing  in  the  former  is  the  usual  differential  quantity, 

a^,  defined  by  Eq.  (1)  or,  equivalently,  by  “  dAG*/dAG°,  whereas  that 

occurring  in  the  latter  is  an  integral  quantity  defined  by  0^-“  (AG*  -  AG*nt)/AG 

Generally  a^  will  differ  from  a^  when  the  Tafel  plots  are  curved,  i.e.  when 

3 

a  depends  upon  AGC .  This  distinction  between  the  forms  of  Eqs.  (12a)  and  (A10) 
arises  since  the  former  can  be  obtained  by  taking  the  temperature  derivative  of 
the  latter  at  a  constant  Galvani  potential  (i.e.  nonisothermal  cell  potential), 
whereupon  da./dT  4  0. 
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TABLE  I  Dependence  of  Experimental  and  Calculated  Transfer  Coefficients  for 

3+ /  2+ 

Cr(OH»)  on  Anodic  as  compared  to  Cathodic  Overpotential 


Anodic 


±(E-Ef)‘ 


c  d 


Cathodic 


c  e 


Anodic  or  cathodic  overpotential. 

Anodic  experimental  transfer  coefficient  corrected  for  work  terms,  at 
listed  overpotentials;  data  from  ref.  1. 

<» 

'Anodic  transfer  coefficient,  calculated  using  individual  force  constants 
for  Cr11  -0H2  and  Cr1  -0H2  bonds  as  described  in  text. 

1 

Cathodic  experiment  transfer  coefficient  corrected  for  work  terms,  at 
listed  overpotentials;  data  from  ref.  3. 

2 

Cathodic  transfer  coefficient,  calculated  using  individual  force  constants 
for  CrI^^-0H2  and  Cr1I-OH2  bonds  as  described  in  text. 


•'Anodic  or  cathodic  transfer  coefficient  calculated  assuming  equal  force 
constants  in  oxidized  and  reduced  forms  [Eq.  (2)]. 


ASLE  II  Comparison  of  Experimental  and  Calculated  Activation  Parameters  for 
Cr (0^) at  Cathodic  and  Anodic  Overpotentials 


ni  „  a 
-E, 


AH* 

corr 


AH*  . 
calc 


kj  mol 


Cr(OH„),  reduction 
2  0 


Cr^H^^  oxidation 


AS* 

corr 


AS*  1 
calc 


J.  K-1  mol"1 


-350 

63 

64 

68 

80 

-300 

68 

66 

68 

80 

-250 

73 

70 

68 

82 

-200 

78 

74 

66 

85 

-150 

84 

77 

65 

88 

350 

11 

17 

-85 

-90 

450 

9 

15 

t90 

-85 

550 

7 

13 

-92 

-80 

650 

5 

12 

-95 

-78 

750 

4 

10 

-92 

-75 

850 

3 

9 

-90 

-70 

900 

2 

8 

-85 

-66 

SwSWmv 


ootnotes  to  Table  II 


'Overpotential,  where  Enl  is  nonisothermal  cell  electrode  potential  with 
reference  electrode  held  at  25eC,  and  Ef  is  the  formal  potential  for  Cr(OH„)^ 
(-660  mV)  at  25cC.  1 

3+ 

Experimental  work-corrected  "ideal"  enthalpy  of  activation  for  Cr ( OH^ ) ^ 
reduction  (cathodic  overpotentials)  or  CrtO^)^  +  oxidation  (anodic  ^ 

overpotentials).  Former  evaluated  using  1  M  NaClO^  supporting  electrolyte, 
and  latter  using  0.4  M  KPF^. 

Corresponding  calculated  "ideal"  enthalpy  of  activation,  obtained  as  described 
in  text. 

^Work-corrected  "ideal"  entropy  of  activation  for  CxiOR^)^  reduction  (cathodic 

overpotentials)  or  Cx(,0R^)^+  oxidation  (anodic  overpotentials),  obtained 

from  corresponding  values  of  AH*  and  k  as  outlined  in  text. 

corr  corr 

Corresponding  calculated  "ideal"  entropy  of  activation,  obtained  as  outlined 
in  text . 


3+ 


Plots  of  the  logarithm  of  the  rate  constant,  log  k,  for  Cr(OH^), 

2  6 

2+ 

reduction  and  CrCOH^)^  oxidation  at  cathodic  and  anodic  overpotentials , 
respectively,  against  the  electrode  potential,  E.  Solid  curves  are  the 
experimental  work-corrected  Tafel  plots  at  mercury-aqueous  interface,  taken 
from  refs.  1  and  3.  Dotted  line  is  anodic  Tafel  plot  corresponding  to  a 
transfer  coefficient  of  0.50.  (Experimental  cathodic  Tafel  plot  has 
slope  of  0.50.)  Dashed  lines  are  anodic  and  cathodic  Tafel  plots  calculated 
using  individual  force  constants  for  Cr*11-©^  and  Cr^-O^  bonds  (see 
text  for  details) .  Dotted-dashed  curves  are  corresponding  calculated  Tafel 
plots,  but  assuming  that  the  outer-shell  reorganization  energy  equals  zero. 

Figure  2 

As  Fig.  1,  but  for  Eu^^)^"*"^"1". 


Figure  3 

Schematic  representation  of  forward  and  reverse  reorganization  energies, 
and  X  ,  for  case  when  the  force  constants  for  the  product  bonds  are  greater 
than  for  the  reactant  bonds.  The  nuclear  reaction  coordinate  is  given  by 
the  dimensionless  parameter  X  in  Eq.  (5). 

Figure  4 

Calculated  free  energy,  enthalpic,  and  entropic  barriers;  AG*V  AH*,  and 

2+ 

T!S*,  respectively;  for  electrooxidation  as  function  of  free-energy 

driving  force,  AG°  [=  -F(E-Ef) ],  using  Eqs.  (5),  (6),  (8),  (11),  and  (13) 


